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We study the effect of turnover of cross linkers, motors and filaments on the generation of a 
contractile stress in a network of filaments connected by passive crosslinkers and subjected to the 
forces exerted by molecular motors. We perform numerical simulations where filaments are treated 
as rigid rods and molecular motors move fast compared to the timescale of exchange of crosslinkers. 
We show that molecular motors create a contractile stress above a critical number of crosslinkers. 
When passive crosslinkers are allowed to turn over, the stress exerted by the network vanishes, due 
to the formation of clusters. When both filaments and passive crosslinkers turn over, clustering is 
prevented and the network reaches a dynamic contractile steady-state. A maximum stress is reached 
for an optimum ratio of the filament and crosslinker turnover rates. 
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The cell cortical cytoskeleton is essential in processes involving cell shape changes Qi- In the cortex, myosin 
molecular motors are assembled in bipolar filamentous structure which bind to actin filaments and generate forces 
by consuming the chemical energy of the hydrolysis of adenosine triphosphate (ATP). The action of myosin motors 
result in the generation of an activ^ contractile stress, whose spatial distribution in the cortex plays a key role in 
cellular morphogenetic processes [1, Q . 

In liviM cells, passive, active crosslinkers and actin filaments are continuously exchanged between the cortex and the 
cytosol [I|. As a result, cytoskeletal networks can release elastic stresses stored in the network and undergo large-scale 
flows. Significant progress has been achieved trough in vitro studies and theoretical analysis of actomyosin networks 
to understand stress generation in networks with permanent filaments and fixed or unbinding crosslinkers [uni. It 
is unclear however what is the role of turnover in stress generation and how filament networks can simultaneously 
rearrange and exert a permanent internal active stress. In vivo experiments suggest that the rate of turn-over is a 
major determinant of force generation by actomyosin networks [121 Il3|. 

We ask here how the rate of turn-over of passive crosslinkers and actin filaments influence the active stress generated 
by motors in the network. We study a simplified mechanical model for a cytoskeletal network in two dimensions whose 
constituents are turning over (Fig. [l}i). Filaments (actin filaments) and motors (myosin motors) are treated as rigid 
rods. Filaments are assumed to have a polarity represented by the arrows as shown in Fig. [TJ Passive crosslinkers 
are assumed to be point-like and constrain the position of the filaments on which they are attached. Filaments are 
able to rotate freely around the cross linker position and around motor heads. Motor heads exert an active force fm 
with a constant magnitude \fm\ = fo on filaments, oriented toward the reverse direction of the arrow (the minus end 
of actin filaments). 



FIG. 1. (a) Schematic illustration of a 2D cytoskeletal network (left). Filaments are represented by black arrows, motors by 
red bars and cross linkers by green dots. The network exerts a stress a on the boundaries of the box, arising from tensions 
acting within the filaments (f/) and within the motors (fm). Motors move towards the arrowhead of filaments (for actin, the 
arrowhead corresponds to the barbed end). Network components stochastically exchange with a reservoir (right), (b) Motors 
move on filaments and rearrange the network on a timescale r. Crosslinkers and actin filaments turn over on longer timescales 
Tc > T and Tf > T. 

To obtain forces acting on filaments, we introduce the effective mechanical potential, U = W+X-g{xf i, ^m,k^ 'iT'm.k); 
where W = — fo %) is the work due to motor active forces, with Ski the position of the fe-th motor head on the 
i-th filament relative to the filament centre of mass, and the sum i) performed for all the pairs of filaments 
(j) and motors (k) connected with each other. In addition, geometrical constraints arise from the conditions that 
cross linkers and motor heads are firmly attached to filaments, and that motor filaments have a fixed length. The 
coefficients A are Lagrange multipliers imposing these constraints, g{xf,i,nf^i,Xm,k,'Hm,k) = O7 where xj^i and nf^i 
(resp. Xm,k, 'i^m,k) indicate the centre of mass position and orientation unit vector of the i-th filament (resp. fc-th 
motor). The motion of motors is taken into account by writing that the position of attachment of the /c-th motor Ski 
relative to the centre of mass of the filament i follows the dynamic equation 

dski _ dU 
dt dski 


( 1 ) 
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with /i a scalar friction coefhcient arising from translational friction between the motor heads and the filament . 

Motors have a typical velocity Vm = and we introduce a reference timescale t = We neglect viscous 

forces arising from the fluid around the network compared to motor-filament friction, and the position and orientation 
of hlaments nf^i are relaxed instantaneously. 

To fix ideas, cortical actomyosin networks in a cell have a mesh size ^ ~ 20 — 250nm [l3, Si and the typical cell 
diameter is several tens of micrometers. We therefore expect actin filaments to have a length If of order 0.1 — l/rm, 
smaller than their persistence length Ip = 16^m 0- Myosins move on actin filaments with velocity Vm ~ 0.1 — S/un/s 
[S, [2l|. The characteristic time for the myosins to move on a filament is r ^ If /vm ^ 0.03 — 10s. To compare this 
timescale to the effect of viscous stresses arising in the solvent of viscosity rj, we note that the velocity of a filament in 
the network subjected to a force /o is ^ /o ln(^/r/)/(? 7 //), giving a timescale ^ /j? 7 /(/o ln(^/r/)), with rf ~ 5nm 
the radius of a filament. Taking the force exerted by a myosin /o 6 — 12piV, ry ~ 10“^Pa.s for water, we find 
Ty ^ 10“^s <C T. Finally, we expect the characteristic times Tc and Tf of the crosslinker and actin filament turnover, 
respectively, to be of the order of 10 s- 1 min [22|, [2^, slow compared to these timescales, such that Tc 2> r and 

Tf » T. 

a. Configuration with one motor and two filaments. We start by discussing the dynamics of configurations in¬ 
volving one motor attached to two filaments (Fig. [5^). The two filaments can be connected to a crosslinker, itself 
connected to the external network, which we consider here to be hxed in space. We distinguish several basic possible 
configurations, depending on the number and positions of attached passive crosslinkers (Fig. [2 )d): both filaments 
can be completely free (i), the two filaments can be crosslinked to each other (ii), or they can be crosslinked to the 
external network (iii)-(v). 

To investigate how molecular motors modify the filaments organisation, we study the relaxation to final state of 
these different configurations. By averaging over possible initial configurations, we evaluate whether motors form on 
average positive or negative force dipoles in the network (Fig. [2]:). The motor force dipole is d = —(/m/2)n- (fi — £ 2 ), 
with n the unit vector giving the motor orientation, and fi and £2 are the forces by which the motors pull or push 
the filaments (Fig. [2Jl). In case (v), £1 and £2 include not only the forces exerted by the motor themselves but also 
the forces originating from the geometrical constraint (Appendix!^. We find that possible initial configurations can 
be classified in 2 categories (Fig. [2]3). In cases (i) and (iv), filaments are moved relative to each other until the motor 
detaches, so that a motor-induced force dipole acts on a transient time t. This force dipole is expansile on average. In 
the quasi-static limit where r <g; Tc, the contribution of transient filament-motor configuration to the overall stress is 
negligible compared to steady motor configurations. When filaments either (ii) have only one fixed attachment point 
to the external network, (hi) are connected to each other by a cross linker, or (v) have both more than 2 cross linkers, 
they relax to a steady configuration where the motor exerts a constant force dipole. By averaging the resulting force 
dipole over possible initial configurations of the two filaments and the motor (Appendix]^, we find that the motor 
exerts a contractile force dipole on average in cases (iii) and (v) (Fig. [2}:). As in Ref. [3], the bias towards contractile 
states arises from instabilities of expansile configurations (Appendix]^. The average force dipole {d) vanishes for 
point-like motors, Im ^ ^ [ 13 ] ■ 

b. Numerical simulations of networks with turnover. We next numerically simulated a network of Nf filaments, 

Nm motor rods and W passive crosslinkers in a square box of width W with periodic boundaries (Fig. [1}. The 
frame of the box is not allowed to deform. We fixed the normalised motor density l^c to 1, where c = Nm/W^ is the 
motor density. To make numerical simulations easier, the geometrical constraints g were replaced by linear springs 
mimicking the contacts at the junctions between motors and hlaments, and at crosslinking points. Initial conditions 
are obtained by randomly positioning hlaments in the network, and timescales are normalized to the reference time 
T. Turnover is introduced by stochastically removing cross linkers, motors and hlaments from the network with 
rates and 1/t^. For simplicity, turnover rates of crosslinkers are taken here independent of the forces they 

sustain [n|. Filaments, motors and passive crosslinkers are added in the network with on-rates fc™ and New 
hlaments take random positions and random orientation, motors take a randomly chosen position on two hlaments 
points separated by a distance Im, and passive crosslinkers are put on a randomly chosen hlament intersection. 

We evaluate the components of the stress tensor (Ty {i,j = x or y) acting on the boundary of the box (Fig. [TJi). 

f f 

The total stress is given by two contributions cjy- = crF -|- a"-, with crF obtained by summing forces acting both 
within the hlaments and ct™ from forces acting within motor rods crossing the boundary of the box (Fig. [1}. In a 
linear elastic or viscous material, crh = 0 in the absence of large-scale deformation. In a non-linearly elastic material 
however, a{j ^ 0 (Appendix [B|) . At steady state, the average stress acting within the motors is given by (Ref. [ 13 ] 
and Appedix la 

cr'” ~ {dmujc) , (2) 

where n, d and c denote the orientation, motor-induced force dipole strength and concentration of bound motors. 
When the dipoles are isotropically oriented, cr™ = with ~ {dc)l2. 
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FIG. 2. Contractile and expansile configurations for one motor and two filaments, (a) Two filaments are attached by fixed cross 
linkers to the external network, (b) (Left) Possible initial configurations for one motor, two filaments and one or several cross 
linkers, and final configurations, (c) Average strength of the contractile force dipole generated by the molecular motor (d) in 
the final configurations (ii), (iii), (v) where the motor does not detach, as a function of the filament length If. (d) Contractile 
(force dipole d > 0) and expansile (d < 0) configurations. 


We first performed simulations without turnover of filaments. We focus on the isotropic component of the stress 
a = {(7xx + o’yy)/2. Figure |3K (blue curve) shows the typical time evolution of the stress, for permanent crosslinkers 
and for crosslinkers with a finite lifetime. Without crosslinker turnover, the system reaches a stationary state where 
the stress a fluctuates around a finite positive value. Figure [5)3 shows the average value of the resulting stress 
(circles) as a function of the number of cross linkers W- For a small number of crosslinkers (W —t 0), no stress 
is observed, tr = 0. Above a critical value of the number of crosslinkers W > N* ~ Nf, a transition occurs and 
a positive contractile stress appears in the system. When W ~ N* but W < N*, although a portion of motors 
generate contractile force dipoles, filaments aggregate in clusters and hence cr ^ 0, suggesting that the transition is 
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associated to the network connectivity (Supp movies Ml and M2). The stress then further increases with the number 
of crosslinkers and eventually saturates to a positive value. Such a transition to contractility as a function of the 
number of crosslinkers has been reported in in vitro reconstituted networks [ifli . , as well as network clustering . 

The average stress only within the motors cr™ is also plotted in Fig. [HJr (cross marks). A similar transition occurs 
for a critical value of the number of cross linkers. The stress ct"* is however larger than the total stress for a large 
number of cross linker Nc > Nf, implying that the filament network is under compression {a^ < 0). For large Nc, 
conhgurations (v) dominate in the network (Fig. [31:). The saturating value of (T™/cro 0.6 is nevertheless smaller 
than the average force dipole obtained by averaging all possible filament orientations (d/2)/(/o/m) 1-0 (Fig. |2l:); 

this is because some configurations do not relax to equilibrium on the characteristic timescale of myosin turnover. 

When crosslinkers are allowed to turn over, the average stress first reaches a positive value before decaying to zero 
(Fig. |31 l, red and green curves), even though the stress within the motors ct"* is still non zero (Fig. |511). The decay of 
the stress correlates with the collapse of the network in an isolated cluster, where filaments accumulate (Fig. jS!;, Supp 
movie M3). To evaluate the decay timescale Tgt, we fitted an exponentially decreasing function a(t) = Aexp(—t/rgt) 
to the simulation results with fitting parameters A and Tst ■ The decay timescale of the stress increases exponentially 
with Nc for large W (Fig. |3f). The relaxation timescale of the stress can be understood as follows: the timescale Tst 
corresponds to the relaxation Maxwell time on which the network becomes fluid, due to turnover of passive crosslinkers 
enabling network rearrangements. This time can be estimated by 



*-l) 


(3) 


with n* = 2Nc/Nf. To obtain Eq. Q, we assume that filaments with at least one crosslinker are fixed, and only 
filaments with no crosslinker attachment can rearrange. We then compute the first passage time at which a filament 
is free from cross linkers, starting from a configuration where the filament is attached with only one crosslinker 
(Appendix |C]). Equation (I3|) accounts for the characteristic time of stress decay in the network for large values of 
Nc/Nf (Fig. H). 

We now turn to simulations where both filaments and crosslinkers turn over. In the cell, actin filaments polymerise 
and depolymerise. Here, we account for this process by simply introducing a rate of filament turnover, tJ^ . Passive 
crosslinkers and motor heads are removed together with the filaments to which they are attached. Remarkably, with 
both crosslinker and filament turnover, the network evolves towards a steady state with a non-zero positive stress 
(Fig. H^-b). As in the previous case, a transition from a non-contractile to a contractile network appears when the 
number of crosslinkers Nc is increased (Fig. HJ:, Supp movies M4 and M5). Note that the stress tr deviates again 
from the average stress within the motors tr™ (Fig. [4}:, inset). The dependency of the fraction of configurations as a 
function of W is qualitatively similar to the case without filament turnover (Fig. HJl). 

We then varied the filament turnover timescale r/ (Fig. |ll:,e,f). We find that the stress reaches a maximum for 
intermediate values of the filament turnover timescale (broken line in Fig. [4^, and Fig. |4f): for slow filament turnover 
Tf ^ Tst, the network collapses in clusters due to crosslinker turnover, while for fast filament turnover Tf Tgt, 
filaments are removed before motors reach a configuration where they can exert a force dipole in the network. We 
find that the optimum value of the ratio of turnover time scale in Fig. 4f, t//tc ~ 10, is of the order of the ratio 
experimentally measured values of turn-over of actin and crosslinker in the cell cortex, rj ~ 15 — 45s and Tc ~ 7 — 14s 

i- 

c. Conclusion. We propose the following dynamic picture for the stress generated in a rearranging network: 
molecular motors move on filaments on a timescale r, fast compared to the crosslinker turnover Tc- When molecular 
motors bind to pairs of filaments, they either displace them and detach, or find steady configurations where they 
generate predominantly contractile force dipoles (Figs. |3| and |5^-c). In networks with permanent crosslinkers, a 
transition to contractile state of the network occurs for a large enough number of crosslinkers (Fig. |3 Jd). When 
crosslinkers are allowed to turnover, the network can rearrange and flow, filaments collapse in a cluster, and the 
total stress in the network vanishes even though the stress only within the motors is still contractile (Fig. |3^,d-f). 
Introducing filament turnover occurring on a timescale r/ comparable to the crosslinker turnover timescale Tc prevents 
this clustering mechanism, and allows the network to generate a steady-state contractile stress (Fig. |3|). It will be 
interesting to investigate stress generation in in vitro experiments where crosslinkers and filaments are allowed to turn 
over. 
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FIG. 3. Stress in a network of motors and filaments, with and without crosslinker turnover, (a) Time evolution of the 
normalised stress cr/ao for = 0, Nc = 1.2Nf (blue), Tc = lOOr, Nc = l.GNf (red) and Tc = lOOr, Nc = 1.2A^/ (green). 
((To = folwiNm/W^) (b) Steady-state isotropic stress a (circles) and motor stress a"' (cross marks) as a function of the 
crosslinkers number Nc for = 0. (c) Fraction of different configurations (i)-(v) in Fig. [2] at steady state as a function of 
the crosslinker number Nc for =0. (d) Isotropic stress (t as a function of the crosslinker number Nc, for finite crosslinker 
lifetime Tc — lOOr, and for several simulation times (rsim. = 300,600,1,200 and 2,400r). Stress within the motors am is 
shown in the inset for Taim. = 300 and 2, 400r. (e) Snapshots of a simulated network with crosslinker turnover {Nc = 1.2Nf, 
Tc = lOOr). (f) Stress decay time Tat as a function of the crosslinker number Nc with crosslinker turnover (tc = lOOr). Solid 
line, theoretical prediction (see main text). Other parameters: Nf = 1,000, Nm = 100, tJ^ = 0, k^^Tc = 20, If = 2lm, 
W = 101m, Tm = lOOr and = 20. 


Appendix A: Calculation of the average force dipole for two filaments 


We detail here how we obtain the average force dipole exerted by a motor on two hlaments in the conhgurations 
shown in Fig. 2. As mentioned in the text, we assume here that crosslinkers are rigidly fixed to the external network. 
The geometry of the two filaments and motor can be described the position of the motors on the two filaments Si, 
S 2 and by the two angles ■i/'i between the motor filament and f-th filament (i = 1,2). The motor positions Si evolve 







































































FIG. 4. Stress in a network with filament turnover, (a) Time evolution of the stress in a network with crosslinkers and filament 
turning-over for r/ = lOOr, Tc = lOOr and Nc/Nf — 1.6. (b) Snapshots of the network evolution corresponding to (a), (c) 
Steady-state stress cr as a function of the number of crosslinkers for Tc = Tf = lOOr (circles), Tc = Tf = lOr (triangles) and 
Tc = Tf = It (squares). Inset: steady-state stress a (circles) and motor stress cr"* (cross marks) for Tc = r/ = lOOr. (d) Fraction 
of different configurations (i)-(v) at steady state, as a function of Nc {tc = lOOr, Tf = lOOr). (e) Heat map for contractile 
stress (T as a function of Nc. and Tc/t/ {tc = t). The broken line indicates Tf = Tst. (f) Stress as a function of the ratio of 
filament and crosslinker turnover rate (r^ = t, Nc = 1.6Nf). Other parameters as in Fig. [3] except = 0. 


according to 


dsi dU 

^ dt dsi ’ 


while the angles Tpi relax instantaneously, so that dU/dipi = 0. 


(Al) 


1. No crosslinkers attached, or one filament is attached by only one crosslinker and the other filament is 

rigidly attached to the external network 


Case (i) and (iv): In these situations, no steady-state is reached as the motor always unbinds from the two filaments. 











































































FIG. 5. Final configurations reached by two filaments bound by a motor, with the filaments attached to a rigid external network 
with a varying number of cross linkers. A, D, G, example configurations for each case, according to the number of cross linkers; 
B, E, H, corresponding possible steady-state configurations; G,F,I, Phase plot of the dynamics of the motor according to the 
coordinates of its position on the two filaments, xi, X 2 - Red lines and dots correspond to stable steady-state configurations, 
while dotted black fines correspond to unstable steady-state configurations. Green regions represent values of the coordinates 
x\ and X 2 which are not geometrically accessible. 


2. Two filaments attached to each other by a crosslinker 

Case (ii): When the two filaments are attached by a single crosslinker, the mechanical potential U reads 


U = -fo{xi -I- 0 : 2 ) -b A 



2xiX2 cos 9 — I 



(A2) 


with xi and X 2 the distances between the myosin attachment points to the crossing point of the two filaments, taken 
positive in the direction of the filament towards which motors move (Fig. [SJA). 6 is the angle between the two filaments, 
such that ni.n 2 = cos 9, with the unit vector giving the orientation of the filament i. In addition, A is a Lagrange 
multiplier ensuring that the length Im of the motor is fixed. Three situations can then occur: 

• when 9 0 and 9 ^ tt and the filaments are not parallel neither antiparallel, the angle 9 is given by 


cos 9 = 


-I- 7-2 — /2 
Xi i- X2 Irn 

2X1X2 


(A3) 


and is free to adjust as the motor position changes. The evolution of the position of the motor on the filaments 
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is given by 


M^=/o (A4) 

^~dt ^ 


such that the motor moves until it detaches from the hlaments, or until the filaments become parallel or 
antiparallel. 

• When the two filaments are parallel and point in the same direction, 0 = 0 and {xi — a; 2 )^ = The center of 
mass of the motor at position (xi + X 2 )/2 follows the dynamic equation 


d Xi+ X 2 
dt 2 


= /o 


(A6) 


and the motor moves on the filament until it detaches. 


• When the two filaments are antiparallel, 9 = tt and (xi + x 2 )'^ = Two configurations are then possible, 
according to whether the two filaments point away or towards the center of the motor from their attachment 
point to the motor. These two configurations correspond respectively to xi +X 2 = —Im (expansile configuration) 
and xi + X 2 = Im (contractile configuration). Any position of the motor on the filament is then a steady-state 
solution, as long as the two motor ends each bind on the filaments. 

To test for the stability of these two configurations, we consider a slight change of the angle between the two 
filaments, 9 = ir + 69, with |(50| <C 1. We take the initial position of the motor to be xi = xq + and 

X 2 = —xo -l- xl‘rn/‘2., with X = for an expansile configuration and y = 1 for a contractile configuration. 
xq = {xi — X 2 )l^ is the position of the center of mass of the motor, relative to the crosslinker joining the two 
filaments, measured positively in the direction of the first filament. We find then the following dynamic equation 
for 69'. 




d69 dx\ 89 




(A7) 


where we have used Eqs. (IA3I) . (IA4I) and (IA5I) from the first to the second line. The sign of the right hand side 
ofEq. dSl) indicates whether the angle 9 increases or decreases when the two filaments are slightly rotated 
away from their antiparallel configuration. Therefore, the associated configuration is unstable when the sign is 
positive, and stable otherwise. We find therefore that the stability of the configuration depends on the position 
of the center of mass of the motor: for jxol < ^m/2 (motor near the crosslinker), the expansile configuration 
is stable, and the contractile configuration unstable. For |xo| > /m/2 (motor away from the cross-linker), the 
expansile conhguration is unstable and the contractile conhguration is stable. 

The flow diagram in Fig. [5]C case (ii) shows the corresponding dynamics in the space of the motor position (xi,X 2 )- 
Red segments indicate the stable steady states. In the green regions, no value of the angle 9 allows for the motor to 
have position (xi, X 2 ) on the two filaments. 


3. Two filaments, each attached by a crosslinker to the external network 

The mechanical potential reads in that case 


U = -fo{xi -I- 0 : 2 ) -I- A 


^Jxl+ xl+l^+ 2xilm COS^i -I- 2x2lm COS ^>2 + 2xiX2 COs('!/'l -|- V' 2 ) — « 


(AS) 


with xi and X 2 the distances between the myosin attachment point and the crosslinker position on each filament, 
and a the distance between the two cross linkers (Fig. 03). A is a Lagrange multiplier imposing that the length of 
the motor is equal to Im- In this subsection, we use for convenience the dynamics of the angle between motor and 
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filaments ipi and ^/’2 to characterise the orientation of the filaments. The dynamics in the limit where ipi and ip 2 relax 
quasi statically is given by 


dx A 

= fo + - {Xi + Im COS ill + X2 COS('01 +^ 2 )) 

dx A 

= fo + -iX 2 + Im COS tjj2 + Xl COs('01 +^ 2 )) 

d'tp A 

= -(xi/msin '01 + a;iX 2 sin(' 0 i +' 02 )) 

ot a 

n 5'02 ^ / , . , ■ , , , 

efia —— =-(x2/msm'02 + a;iX2 sm(' 0 i +02)) 

ot a 


(A9) 

(AlO) 

(All) 

(A12) 


where e <C 1 is a factor that vanishes in the quasi-static limit where the filaments rotate adiabatically. Solving for the 
Lagrange multiplier A by imposing the constraint that the length of the motor is equal to Im, we obtain: 


3(0:1 + a;2)(l + cos(0i+02)) +/m(cos0i + COS02) 

^ - 

A = {xilm sin 01 + a;ia;2 sin(0i + 02))^ + {x 2 lm sin 02 + a;ia;2 sin(0i + 02))^ 

B = {Xi + Im COS 01 + 0:2 COS(01 + 02))^ + (0:2 + Im COS 02 + 0:1 COs(0i + 02))^ 


(A13) 


where the coefficient A vanishes when the filaments are aligned with the motor. Several situations can again be 
distinguished: 

• When the filaments are neither parallel nor antiparallel, from Eq. (IA13(1 . A ^ e in the limit e —>■ 0. Eqs. (IA9I1 
and (|A10I) then yield at the lowest order in e 


= /o + 0{e) (A14) 

= /o + 0{e) (A15) 

such that no steady-state exists in that situation. 

• When the filaments are parallel and point in the same direction, (0i,02) = (0,7r) or (0i,02) = ('7r,0). In that 
case A = 0 and no steady-state exists for the motor, which runs on the two filaments before detachment. 

• When the filaments are antiparallel, (0i,02) = (0,0) (expansile configuration) or (0i,02) = ('7r,7r) (contractile 
configuration). In the expansile configuration, xi + X 2 +lm = with y = ±1. In the contractile configuration, 
Xl + X 2 — Im = XO', with the same rule applying for x- In both cases, Eq. (IA13I) yields A = —x/o and any 
position of the motor is a steady-state. 

To test for the stability of these steady states, we consider the dynamics of 0i and 02 around 0i = 02 = 0* 
with 0* = 0 or TT, 01 = 0* + (501 and 02 = 0* + 502- To regularise the dynamics around the parallel filaments 
state, we consider a situation with finite e and consider perturbations verifying (50 -C \/e, such that A ~ — x/o- 
From Eqs. (jAllIl and (IA12I) . we have then 


efj,a [{xilm COS0* + XiX2)S'ijji + X1X2S1IJ2] 


e^a- 


dt a 

, dS'ip2 ^ ^ 

dt a 


[xiX2Stpi + {X2lm COS0* + Xia;2)(502] 


(A16) 

(A17) 


so that the state 0i = 02 = 0* is stable if xixilmcosijj* + X 2 lmCOStp* + 2 x 10 : 2 ) < 0 and {xilmcostp* + 
xiX 2 )(x 2 Zm COS0* + X 1 X 2 ) — X 1 X 2 > 0, whereas it is unstable otherwise. Taking the initial condition to be 
Xl = xo/2 + (xa — Im cos 0*)/2 and X 2 = — xo/2 + (xa — Im cos 0*)/2 at the steady state, the stability condition 
is satisfied when 


x(xq +lm~ «^) > 0 and X1X2XCOS0* > 0 . (A-18) 

The results are summarized in the flow diagram in Fig. [Sp. Red segments indicate the stable steady states. In the 
green region, there are no possible values of 0i and 02 for given values of xi and X2. 
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4. Two filaments rigidly attached to the external network 

This situation corresponds to case (v). The mechanical potential U reads as for case (ii) 


U = -fo{xi + 2 : 2 ) + A (y\lx\+xl- 2 xiX2 cos 9 - , 


(A19) 


with xi and X 2 the distances to the crossing point of the two hlaments, taken positive in the direction of the filament 
towards which motors move (Fig. [5)2). 9 is the angle between the two hlaments, such that ni.n 2 = cos0, with 
the unit vector giving the orientation of the hlament i. A is a Lagrange multiplier ensuring that the length of the 
motor is hxed. The dynamic equation for the motion of the motor on the two hlaments Eq. (EH) then reads 


9'-^ = fo - 7 —(a;i - X 2 cos9) 


dt 

dx2 


In 

A 


H— = fo - “( 2^2 - xi cos 6 »), 


(A20) 

(A21) 


where d is a hxed angle, and the Lagrange multiplier A is obtained from the constraint that x^+x^ — 2xiX2 cos 9 = l^: 

y folm{xi+X2)il-COs9) 

{xi — X 2 cos 0 )^ + (x 2 — a;i cos 9)^ 

Two solutions can be found for the motor position on the two hlament, assuming that they are not parallel (0 7 ^ 0 
and 9 tt): 


Xi = X2 = ± 


2 sin) 


(A23) 


As pointed out in Ref. Q , a linear stability analysis around the steady-state indicates that only the positive solutions 
xi = X 2 = /rn/[2sin(0/2)] is stable. The corresponding force dipole at equilibrium is positive and is given by 

^m/o 


d = 


Sind)- 

Therefore, the dipole formed on two rigidly hxed, non-parallel hlaments is always contractile (Fig. [SH-I)- 


(A24) 


5. Averaging force dipole 


To obtain the associated average force dipole exerted by the motor in these different conhgurations, we proceed as 
follows: denoting as 0 ^ 2 ^ the initial distance between the two hlaments, 0 ^*] and 0^*2 the two angles giving the initial 

orientations of the hlaments, x^^^, x^p 2 initial position of the two motors on the hlaments, and Xk (/c = 1 , • • • , Nc 
with Nc the number of crosslinkers attached to the external network) the position of the crosslinkers on the hlaments, 
we compute the motor-induced force dipole d reached at steady-state as shown in Fig. [5] cases (ii), (iii) and (v) by 


d{c^°^ ,9fl,9't^,x 


1 ( 0 ) 

/, 2 ’ 




^mfo 


^mfo 


2 COS 'l/ji 2 COS 2^2 


(A25) 


where the angles ipi and V '2 between the motor and two hlaments at steady-state are functions of cJ 2 ^ 
^ 2 °^ and {rfcjfe. Equation (IA25I) is obtained from the dehnition of the force dipole 


d=-^n.(fi-f2) 


(A26) 


with n the unit vector giving the motor orientation pointing towards hlament 1 , and with i = 1,2 are the forces 
exerted by the motor on the two hlaments. The dependencies 1/ costpi (i = 1, 2) in Eq. ()A25I) arise from the condition 
that the projections of the forces on the hlament i have magnitude fo (Fig. (B]). 

The average force dipole is then obtained by 


< d>= 




d{x 


(0) 

1 ! 




(0) fl(0) fl(0) ^(0) 


-d(N^> 9^^> 9 

f,2- 


X 




{rk}k) , 


(A27) 


where the integration J x^^^) runs for all possible initial motor conhgurations x^°\ x^^ that eventually reach 

a steady state, and F = J (i(xj°\ X 2 °^)l. 
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FIG. 6. Schematic of the force dipole exerted by a motor. Each motor exerts a force of magnitude /o parallel to the filament. The 
total force exerted by the motor on the filament are equal and opposite at steady-state and denoted /i and / 2 . The contribution 
of the force normal to the filament arises from geometrical constrains. Overall, this results in a tension / = |fi| = |f 2 | acting 
within the motor, which is used here to define the force dipole exerted by a motor. 


Appendix B: Relationship between the stress and the 


average force dipole 


In this section, we discuss the two contributions to the stress generated by the filament and motor network, cr(j and 

a^. We point out that in the absence of large scale deformations, a(j can contribute to the total stress, when the 
elastic response of the filament network is non-linear. 


1. Two-dimensional, linear elastic material 


We consider here a two-dimensional elastic material subjected to the force dipoles exerted by motor filaments. The 
motor force dipoles consist of two opposite forces = ±/^n^, separated by a distance Im, and where the index k 
label the motors. The motors induce a deformation in the elastic material. The total stress in the system is then the 
sum of the resulting elastic stress, denoted a-j, and the stress generated within the motors, denoted cr^: 

(Bl) 


We start by discussing the average stress created by an ensemble of motors in a 2D material. Following Ref. [^ . 
we now show that when the system is homogeneous, force balance on a section S of the network enclosed in a contour 
C allows to relate the stress generated within motors to the concentration and orientation of motors. To evaluate the 
stress, we first note that each motor k corresponds to a line under tension with length Im, orientation n^, centered 
at position x^, and joining the two points at coordinates — lm/2n^ and x^ + lm/2n^ (Fig. [Zl^). Therefore, the 
two-dimensional stress field within the motors can be written 

im 

= X] - (x'^ -b n'^u)) (B2) 

where u is a coordinate going along the line under tension /^. The average stress within a region S of surface area S 
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B 



|yWv|AAA/^|^AAA/^|^^AA/|AAA/^^ 

< - > 

f f ' 


> X 


Motors are 
turned on 


If 

0'-^|vAAA/N|\AAA/\|vW\/jwWs|^^ 

I 

FIG. 7. A. Schematic of a 2D elastic material, subjected to forces arising from force dipoles exerted by motor filaments. Motors 
result in a stress aij acting at the boundary of the box. B. Schematic of a ID chain of non linear springs subjected to force 
dipoles. The chain contains N springs and is fixed at both ends. Motors result in a force a along the chain. 


is then given by 


^ J^dyia^{x) 

im 

^ ^ [ , du6{x - + n^u)) 


Sjs 




(B3) 

(B4) 

(B5) 


= c{dninj) (B6) 

where c is the concentration of motors, df = Imf^ is the dipole strength of motor fc, and the averaging (•) is performed 
over space. 

In a linearly elastic material, the average stress within the network of hlaments is given by 


(T^j — 2E KUkkdij 


(B7) 


with E and K a shear and bulk elastic moduli, and tty = ^{diUj + djUi) the gradient of deformation. The average 
stress generated in the elastic material in a region S with contour C and surface area S is then given by 




2E f , K-E 


dxUkkdii 


S 

E f n E f K-E f 

= — y diViUj + -g J dlUjUi H- - — J dlVkUkd 




(B8) 

(B9) 

(BIO) 
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where u is the vector normal to the contour C and dl an infinitesimal line element on the contour. If we consider a 
square box whose boundaries are fixed (Fig. [TjA.), or periodic boundary conditions, the contour integrals in Eq. (IBIOI) 
vanish. Therefore, for a linear elastic material with fixed boundaries, the average stress arises entirely from the forces 
acting within the motors. This however does not apply to a non-linearly elastic material, as we show in the next 
section. 


2. One-dimensional chain of non-linear springs 


We consider a simpler example in ID of a periodic chain of N elastic springs. Each spring is located between 
positions Xi and Xi+i, with initial resting position Xi = il. The two points at the end of the chain, xq and xn, are 
not allowed to move. In addition, Nm motors are acting in parallel to a fraction n = N^/N of the springs (Fig. [Tjl). 
The motor exerts a constant force /. The springs have a non-linear force-extension relation 


fe — k ——h ^2 



(BII) 


with /e the force exerted by the spring, Al = Xi+i — Xi — I the extension of the spring, and k and k 2 are two spring 
constants. We assume that the springs are weakly non-linear, k 2 /k 1. In the initial resting position, Al = 0. The 
motors are then turned-on, driving a deformation of the the springs in the chain. The contraction of the springs 
which are in parallel with the motors is denoted —Aim- Because the overall length of the chain is kept fixed, the 
deformation of the free springs is then given by Almn/{1 — n). We denote by a the total force acting within the chain. 
Force balance imposes that the total force within the springs and the motors is fixed and equal to cr, giving 


, Aim , 

O’ = —k — -h k2 


Aim \ , r 1 ^ 

—) 




/ /\lm 

\ I 1 — n 


2 


(BI2) 


where the second part of the equality is the total force within the springs in parallel with a motor, and the third part 
the force within the free springs. Solving these equations and expanding to first order in fk 2 /k^ <C I, one obtains 


a = fn- 


k 2 pn{l - n) 


or using the concentration of motors c = n/l, 


a = flc H—|^cZ(I — d) 

fC 

= a^ + af 


(BIS) 


(BI4) 

(BI5) 


where cr"* = flc is the average tension exerted within the motors. To leading order in the spring non linearity, k 2 —t 0, 
the force within the chain reduces to the average motor tension cr”*, in agreement with Eq. (IB6I) . The non-linear 
elastic behaviour however brings a correction to this term cr^ = k 2 f^cl{l — cl)/k‘^, proportional to the motor force 
squared, /^. 


Appendix C: Network relaxation time 




Filament 


FIG. 8. Schematic of a filament in the network. Crosslinkers bind to the filament with rate and unbind with rate 1 /tc. 


We derive here an approximate expression for the network characteristic relaxation time. We consider a filament 
within the network, crossing other filaments that are themselves immobilised. We expect this last assumption to 
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be valid for a large enough number of crosslinkers. The number of crossing points is assumed to largely exceed the 
total number of crosslinkers in the network. Crosslinkers bind and unbind the filament at crossing points with other 
filaments. A filament with two attached crosslinkers is completely fixed, and can only rotate when it has one attached 
crosslinker. For simplicity we consider here that the filament can not rearrange to relax stresses unless no crosslinker 
attaches it to other filaments in the network. Once a filament is free, it can move until a crosslinker binds to it and 
immobilise it. We therefore estimate the network relaxation time as the mean first passage time to a state where no 
crosslinkers bind the filament, from a state where one crosslinker binds the filament. 

We consider the probability of having crosslinkers on the filament, P(nc). New crosslinkers bind to the filament 
with rate as crosslinkers bind to all filaments in the network with rate (Fig. [5|), 

and each crosslinker binds two filaments. In addition, crosslinkers unbind from the filament with rate (Fig. [8|). 

The probability P{nc) then follows the master equation: 


= "^Pinc + 1) + fclP(nc - 1) - + Ki ) Pine) 


(Cl) 


d. Stationary distribution. The stationary distribution of Eq. m is a Poisson distribution: 


P{nc) = e 


iKiTcY 


(C2) 


with mean and variance n* = kYnTc — 2Nc/Nf for k^^Tc » 1. 

e. Mean first passage time. We now want to obtain the mean first passage time to reach a state where the 
filament has no bound crosslinker, Uc = 0, starting from a configuration with n crosslinkers attached on the filament. 
We denote T„ this first passage time, and we follow a standard procedure to obtain its value [H, [2^ . T„ satisfies the 
following equation for n > 0: 


= 


k' 


,cf 

on 


Kh + n/r, 


'Tn+l + 


u/Tc 


Uci 

^on 


■ n/r. 


'Tn-l + 


iUcf 

^on 


■ njTc 


(C3) 


where the first term corresponds to the probability of moving to a state with n + 1 bound crosslinkers, times the 
waiting time from the state n + 1, the second term is the product of the probability to move to a state with n — 1 
bound crosslinkers times the waiting time from the state n — 1, and the last term is the average time spent in the 
state n. In addition, Tq = 0 by definition of the first passage time Tn, and we take a reflecting boundary condition at 
infinity, implying Tn — Tn-i —>■ 0 for n —> oo. Equation (IC3I) can be rewritten 


kfniTn - Tn+l) + -{Tn - Tn-l) = 1 


(C4) 


or defining Zn = Tn+i - Tn, 


“ r,cf _ ^"-1 iLcf 

'^on'c ^on 


To solve this equation, one introduces Zn = Zn{kfn'^c )'^which satisfies then: 


— ^n—1 


{KnTcY 


KW. 

Solving Eq. (IC6I) then yields the following expression for and r„, using that Zoo = 0: 

\k 


Z = ^ {^oliT 


k‘Ynk\ 

k=n+l 


— 


{koliTcY 


(^ora'^c)^ 

, ^ Kikl 

k=n+l 


n —1 


Tn = Yl 


kl 


k=o^^oirc)'^ mtY+i 


(C5) 

(C6) 

(C7) 

(C8) 

(C9) 
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From this last expression, one finally obtains Ti, 




(Kl'rcy 


Uci ' 


- 1 ) 



(CIO) 

(Cll) 

(C12) 


The time Ti therefore contains a factor increasing exponentially with the number of crosslinkers per filament n*. 
We use the time Ti as an estimate for the network relaxation time. In Fig. 3f, we verify that the time Ti provides a 
good approximation for the relaxation time of the network. 


Appendix D: Effect of motor concentration and dispersion in motor forces 



FIG. 9. A. Isotropic stress as a function of the number of motors Nm, for two values of myosin turn-over. In these graphs, 
we performed the simulations for the cases with identical motors (red circles) and with the a dispersion of 20% motor force 
and friction (blue triangles). Non-specified parameters are as in Fig. 3 of the main text. B. Isotropic stress as a function of 
the number of motors, in the absence of a cross linker. Green diamond marks, identical motors; no stress is generated. Blue 
triangle marks, a dispersion of 20% motor force and friction results in a non-zero stress for a high enough number of motors. 
Non-specified parameters are as in Fig. 3 of the main text. Stress was averaged here over lOr after initialization. 


We have performed additional simulations measuring the isotropic stress cr as the concentration of myosin motors 
Njn/W^ is varied (Figure. |9]A). The reference stress cto is defined as being proportional to motor concentration, 
such that a ^ ao indicates that the stress is proportional to the myosin concentration. We find that a has a weak 
non-linear dependence as a function of the myosin number Nm for the simulation parameters plotted in Figure [SIA, 
and is nearly linearly increasing with Nm for slow enough turn-over. This indicates that the stress generated by the 
simulated network is roughly proportional to the myosin concentration. 

In our simulations, motors are identical and have the same stall force; as a result no stress is generated in the 
absence of cross linkers (Fig. 3 and Ref [13 )■ Introducing a dispersion in motor friction and stall force however can 
result in stress generation as the number of myosin motors is increased (Figure. [Sjl). The overall stress generated is 
much smaller than stresses generated for the same parameters with a cross linked network. 


Appendix E: Details of the model used for the numerical simulation 
1. Main Components 

We introduce three components in our simulations: actin filaments, crosslinkers, and myosin filaments. 

Myosin filaments are represented by rigid rods with a length Im ■ To represent the attachment between myosin and 
actin filaments, two springs are added at the end of each myosin mini filament with stiffness kms and reference length 
0, such that the tension in the end spring fms is fms = kmslms, with Ims the length of the connecting spring. In 
practice kms is taken to be large enough that the spring maximum extension is very small compared to other lengths 
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in the simulation. The position of a myosin mini filament /c on a filament i, relative to the center of mass of the 
filament i, is denoted s^i- 

Actin filaments are treated as rigid rods with a finite length, If. The position of the centre of mass of the filament 
is denoted Xf^i and the unit vector giving the orientation of filament i is denoted 

Crosslinkers are represented by springs connecting two actin filaments, with a stiffness kx and reference length 0, 
such that the tension with an crosslinker is fx = kxlx with lx the length of the crosslinker. In practice kx is chosen 
large enough such that the length of the cross linker is very small compared to other lengths in the simulation. 

We simulate the network in two dimensions x,y, in a periodic box of size W. Simulations are initialised by 
positioning filaments in the box with random positions and orientations. 


2. Dynamics 

In the simulations, bond myosin minifilaments move on actin filaments. Myosin filaments detach from a filament 
when they reach the end of a filament. In addition, when myosin turnover is taken into account, they can spontaneously 
bind and unbind filaments. Crosslinkers can also bind and unbind filaments when crosslinker turnover is taken into 
account, and filaments can be removed and added in the network when filament turnover is taken into account. At 
every step of myosin motion, the network is relaxed quasi-statically to equilibrium. 


a. Myosin motion 

At every step, the end position of the myosin minifilament k on filament i, denoted Ski, is updated according to the 
following equation: 

dshj , . 

= /o - ffes • Hi. (El) 

where ffcs is the tension of the myosin end-spring connected to the filament, /o is the active force generated by the 
myosin on the filament, and /r is an effective friction coefRcient between the motor and the filament. The equation 
above is discretised with an Euler explicit scheme with time step dt. 


h. Filament and crosslinker turnover 

Crosslinkers are added with a rate k^^ and removed with a rate = I/tc. The number of bound cross linkers 
Neb and unbound cross linkers Ncu add up to the total number of cross linkers = Neb + New At every step, 
each unbound cross linker has a probability of being added to the network k^^dt, and each bound cross linkers has a 
probability to be removed from the network Unbound cross linkers are added by looking randomly for a free 

cross linking point on two filaments, and attaching the crosslinker there. 

Similar rules apply to turnover of myosin and actin filaments. The binding positions of two ends of the myosin 
filament are determined in the following way: Firstly, an actin filament and the position on it are randomly chosen, and 
one end of the myosin filament attached there. Then, the position of the second end of myosin filament is randomly 
chosen from all the attacheable positions, i.e. all the positions on actin filaments located Im away from the position 
on which the other end is attaching. 


c. Quasistatic relaxation 

At every time step, the network is relaxed quasi-statically to equilibrium. Below, we denote t* a fictitious time 
coordinate used for quasi-static relaxation. Quasi-static relaxation is performed by updating the filament center of 
mass Ti and orientation = {cos 6i, sin 9i) according to the following equation: 



k 


r 


dOj 

dt* 


= Gz ■ 



SkUi X 



(E2) 


a- 


(E3) 
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where fik is the force acting on filament i from the crosslinker or motor k. at and are two translational and 
rotational fictitious friction coefficients, used for the quasi-static relaxation, and is the direction orthogonal to the 
plane of simulation. The same equations are used to iterate the position of myosin motors, and orientation n^. 

Iteration is performed by discretising equations (IE2I) and (IE3|) with an Euler explicit scheme, until the system 
reaches quasi-static equilibrium. In practice, a criterion must be used to specify when the system is close enough to 
equilibrium. To do so, the squares of the filament translational and rotational velocities are averaged according to: 



where N is the total number of actin and myosin filaments. The quasi-static relaxation is completed when e is less 
than a threshold parameter Cmax- 


Appendix F: Supplementary movie legends 

• Supp. Movie Ml Simulation of a network with no crosslinker turnover, no filament turnover, Nc = 800 and 

Tm = IOOt. Total simulation time, 1600r. The number of cross linkers is below the threshold for the network 

to exert a contractile stress. 

• Supp. Movie M2 Simulation of a network with no crosslinker turnover, no filament turnover, = 1100 and 

Tm = lOOr. Total simulation time, 1600r. The number of cross linkers is above the threshold for the network 

to exert a contractile stress. 

• Supp. Movie MS Simulation of a network with no filament turnover, crosslinker turnover with Tc = lOOr, 
Nc = 1200 and Tm = lOOr. Total simulation time, 1600r. The network collapses and does not exert a 
contractile stress in steady state. 

• Supp. Movie M4 Simulation of a network with filament turnover with Ta = lOOr, crosslinker turnover with 
Tc = lOOr, Nc = 800 and Tm = lOOr. Total simulation time, IGOOr. The network reaches a steady-state where 
no contractile stress is exerted. 

• Supp. Movie M5 Simulation of a network with filament turnover with = IOOt, crosslinker turnover with 
Tc = IOOt, Nc = 1200 and Tm = IOOt. Total simulation time, 1600t. The network reaches a steady-state where 
a contractile stress is exerted. 
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